NMR evidence for very slow carrier density fluctuations in the organic metal 

(TMTSF) 2 C10 4 
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We have investigated the origin of the large increase in spin-echo decay rates for the 77 Se nu- 
clear spins at temperatures near to T — 30K in the organic superconductor (TMTSF^CICU. The 
measured angular dependence of T 2 -1 demonstrates that the source of the spin-echo decays lies with 
carrier density fluctuations rather than fluctuations in TMTSF molecular orientation. The very long 
time scales are directly associated with the dynamics of the anion ordering occurring at T = 25K, 
and the inhomogeneously broadened spectra at lower temperatures result from finite domain sizes. 
Our results are similar to observations of line-broadening effects associated with charge-ordering 
transitions in quasi-two dimensional organic conductors. 
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(TMTSF) 2 C10 4 is a member of the Bechgaard salts 
family of organic conductors and superconductors. While 
the original discovery of organic superconductivity was 
made in the isostructural system (TMTSF) 2 PF 6 at 
Orsay 1 , the perchlorate salt was interesting because it 
was the only composition undergoing a superconducting 
transition (T c — 1AK) at ambient pressure^. Another 
difference, when compared with the hexaflouride salts is 
that the counterion is non-centrosymmetric. As a result, 
there is a symmetry breaking of the high-temperature 
crystallographic space group (Pmma), when the CIO4 
counterions orientationally order at Tao = 25if with 
wavevector Q = (0, 1/2, 0) 4 . This particular Q leads 
to inequivalent TMTSF stacks and therefore also two 
bands crossing the Fermi energy Ep. Some of the low- 
temperature properties of (TMTSF) 2 C10 4 that differ 
from the PFg salt are attributed to this difference^SiL&. 

The phase transition at Tao is a metal-metal phase 
transition driven by lattice Coulomb effects at which a 
disproportionation, or charge order (CO), of carrier den- 
sity between the adjacent stacks accompanies the coun- 
terion ordering. By now, CO is known to be ubiquitous 
to the analog TMTTF-based salts with counterions rang- 
ing from hexaflourides like AsFg to non-centrosymmetric 
SCN or Re04. Those systems are all insulators, either re- 
sulting from the broken-symmetry of the CO phase tran- 
sition itsel&iSiii or resulting from a dimerization of in- 
trastack bond distances between TMTTF molecules 1 ^* 1 ^. 
In the case of the insulators, there is little doubt that 
electron-electron interactions play a crucial role; never- 
theless there is empirical evidence that coupling to the 
counterions allows the CO transition to take place. We 
examine the effects of the transition on the NMR proper- 
ties as a window into what differences might be observed 
when compared to what is seen in the insulators. 

Here we explore the effects of the anion ordering on the 
77 Se NMR spectrum and relaxation in (TMTSF) 2 C10 4 . 
In this case, the transverse spin relaxation rate is strongly 
affected by slow fluctuations of the electronic carrier den- 
sity above Tao- The slow fluctuations have been ob- 
served previously by Takigawa and Saito (TS) 1 ^, but 
with the difference that our experiments are done us- 
ing single crystals, and we are able to demonstrate that 



the carrier density fluctuations rather than molecular ori- 
cntational fluctuations are responsible. The linewidth 
broadens homogeneously in association with these fluc- 
tuations and narrows on cooling further. From experi- 
ments undertaken below Tao on quenched and relaxed 
samples, we conclude that the spectrum is dominated by 
disorder effects. This is very different from the insula- 
tors, where distinct symmetry-breaking signatures of the 
ordered phase are observed—. Even in the presence of 
the disorder, (TMTSF) 2 C104 is a superconductor, and 
the pairing is probably not s-wave. 

The (TMTSF) 2 C104 crystals were grown using stan- 
dard electrolysis techniques. In this case, the dimen- 
sions are 6.4mm x 1.2mm x 0.8mm and the mass is 
m = 4.2mg. Typically, crystals grow longer along the 
highly-conducting stack direction relative to the other 
directions. Our goal was to investigate the effect on 
the hyperfine fields resulting from the CIO4 anion or- 
dering, so the appropriate coil geometry is unusual when 
compared to most NMR investigations of the Bechgaard 
salts. We expected the hyperfine fields are principally 
dipolar 14 , arising from the p z orbital on the Se atoms, 
extending out of the plane of the molecule. As a result, 
the strongest hyperfine fields are with the dc field aligned 
along the stacking axis, so we chose the symmetry axis 
of the coil to allow for sample rotations with the external 
field Bo lying in the a-c* plane. The AO transition oc- 
curs at Tao = 25K, so we cooled at AOmK/min. through 
it to reach the highly-conducting relaxed state. 

Shown in Fig. [T]is the spin-echo decay rate T 2 ~ (here 
we define T 2 as the time corresponding to 1/e decay) vs. 
temperature for several angles of the applied field B = 
4.91T. As observed by Takigawa and Saito 1 ^, we see 
a substantial increase in the relaxation rate on lowering 
the temperature through 30K . The questions we answer 
are how this increase in T 2 _1 is related to anion ordering 
and the mechanism by which it occurs. By rotating the 
sample, we establish the anisotropy in T^ 1 in the vicinity 
of T = 30K. The size of the peak is found largest for 
the field applied along the stacking axis (Bo || a, 8 = 
0), somewhat smaller for the c*-direction, and even less 
when Bo is rotated away from a by 6 = 6 m = 55° (inset; 
also see Eq. EJbelow, noting that 3cos 2 # m — 1 = 0), where 
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FIG. 1: (Color Online) Temperature dependence of T^ 1 for 
three orientations of the crystal relative to the external mag- 
netic field. The inset illustrates the geometry. 8 indicates the 
angle between the direction of applied field Bo and a axis. 
The open circles are taken with Bo || a {9 = 0), the open 
triangles are taken with Bo || c* (9 = 90°), and the open 
squares are taken with 9 — 8 m . The dotted lines serve as 
guides to the eyes. 



the anisotropic part of the hyperfine field vanishes. At 
the peak, the rate is too fast to allow for a measurement 
with the field along a. 

A more complete angular dependence of the spin-echo 
decay rate is shown in Fig. |2J at T = 38K, where 
the relaxation rate is enhanced over a temperature- 
independent background. We also notice that there is 
a change in linewidth and lineshape on cooling through 
the temperatures where T 2 _1 peaks. In the high tem- 
perature regime (i.e., T rj 50K), the lineshape is very 
close to Gaussian. It changes to Lorentzian on the high- 
temperature side of the peak and returns to Gaussian 
(although considerably broader) on the low-temperature 
side of the peak. 

In what follows, we make the case that the peak in T 2 
vs. T is a result of the slowing down of carrier density 
fluctuations linked to the inequivalence of TMTSF stacks 
developing at the anion ordering transition. The spin- 
echo decay is affected by the density fluctuations through 
the z component of the hyperfine fields, because local 
temporal variations result in dephasing of the precessing 
spins involved in a spin-echo experiment. To be more 
specific, suppose that a nuclear spin is situated in an 
environment in which the local field switches randomly 
between two discrete ±5h z on a time scale t c . If this is 
the only source for the spin-echo decay, then the spin- 
echo amplitude decreases as 

G(t) = G{Q)e-^ 5h ^ T " t , (1) 
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FIG. 2: (Color Online) Angular dependence of T 2 _1 (red open 
circles) at T = 38K. The pink dotted line and the green 
dashed line are the expected angular dependence of T 2 _1 in 
TS model, corresponding to the molecular motion amplitudes 
of 2° and 10°, respectively. 



in the extreme narrowing limit (jSh z T c <C 1), with 7 the 
gyromagnetic ratio. When Eq. Q]applies, the spectrum is 
said to be motionally narrowed. The spin-echo decay rate 
increases when "fSh z r c grows larger, as it would on the 
approach to the anion ordering transition. Dephasing is 
most efficient when "fSh z T c ~ 1, where we expect a peak 
111 IT 1 . Below, we evaluate the form of the hyperfine 
coupling, and use it to describe the angular dependence 
of the spin-echo decay rate. 

The angular dependence of the first moment as well 
as the spin-lattice relaxation rate T 1 ~ 1 at T = 38K are 
shown in Fig. |3| We assume these are consistent with a 
uniaxial hyperfine coupling of the form 

A(9) = A lso + A ax (3cos 2 6-l), (2) 
K(9) = A{0) X S \, (3) 

with v the volume per formula unit and \ s the dimen- 
sionless spin susceptibility. Ai SO is the contribution from 
core polarization of the Se ions. We assume that A ax is 
dominated by the intra-atomic Se p z orbital and takes 
the form A ax — (2/5) < r~ 3 > a with a the carrier den- 
sity on a single ion. The upper limit for a is 0.25; it is 
reasonable to expect less, and TS used a = 0.19 calcu- 
lated in Ref>i&. The stacking (i.e., a) axis is very nearly 
coincident with the principal axis orthogonal to the plane 
of the TMTSF molecule. 

Then the hyperfine field contribution to the spin-lattice 
relaxation rate also varies with angle, even though the 
spin susceptibility is isotropioit For Bq || a {0 = 0) and 
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FIG. 3: (Color Online) First moment of the spectra (solid 
circles) and spin- lattice relaxation rate T^ 1 (solid squares) 
vs. angle 9 at T = 38is:. 



where correlations are ignored, we have 
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with 7 e the gyromagnetic ratio of the free electron. For 
the angular dependence of the relaxation rate, we obtain 
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where a = A iso /A ax . 

Inserting the accepted value for the spin susceptibility 
of x s = 1.6 x 10~ 4 emu/mole(f.u.)H into Eq. |3 together 
with K ax = 10.7 x 1CP 4 extracted from the angular de- 
pendence of the first moment (Fig. [3} , we obtain a value 
for < r~ 3 >= 15<2q 3 , where ao is the Bohr radius. This 
value is larger than the 9.3a, - ' 3 obtained for Se atoms 
from Hartree-Fock calculation*^. The solid line running 
through the relaxation rate data shown in Fig. [jjjis a least 
square fit to Eq. from which we extract a = 0.37 19 . 
Returning to Fig. [2] and using this value of a, we obtain 
the solid line through the data points using the function 



T 2 = p{K{6)~{B q )tc(T) + B 



(6) 



with Bq the external applied magnetic field, and B a tem- 
perature and angle independent constant. The first term 
of Eq. results from local field fluctuations originating 
in variations of the hyperfine coupling (Eq. , which are 
caused by carrier density fluctuations. At higher tem- 
peratures, t c is sufficiently short so as to make hyperfine 
field fluctuations inconsequential for the spin-echo decay, 
and it grows upon cooling with the slowing of the lattice 
fluctuations. (3 is a constant that measures the amplitude 
of the hyperfine field fluctuations relative to the average. 

Although the dependence of the spin-echo decay rate 
on the external field is the same as in TS model, the 
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FIG. 4: (Color Online) Homogeneous and inhomogeneous 
linewidth as a function of temperature, taken with Bo || c*. 
The open circles are the inhomogeneous linewidth of the spec- 
tra, defined as the square root of the second moment. The 
open squares are the homogeneous linewidth, denned as the 
HWHM of the Fourier Transformation of the spin-echo de- 
cay. The triangles are the linewidth taken in the completely 
quenched state. The lines through the data points are guides 
to the eyes. 



dependence on angle is very different to that expected 
for orientational fluctuations of the molecule. For con- 
trast, we show the relative angular dependence of this 
mechanism for orientational fluctuation amplitudes of 
SO = 2°, 10° in Fig. [3 In that case, the spin-echo decay 
rates induced by the local field fluctuations are largest 
where dK/d6 is maximum, namely near to 6 m . Instead, 
we observe a minimum at 9 m . The angular dependence 
is consistent with fluctuations in the hyperfine coupling, 
which includes the factor a representing the fractional 
number of carriers per Se atom. Our interpretation is 
that as the correlation length for the counterion ordering 
grows, so does the characteristic time scale for the fluc- 
tuations in the usual way. In the ordered state there is 
an incquivalent carrier density on adjacent stacks, as a 
result of the period doubling of the lattice in the b direc- 
tion. The fluctuations in the component of the hyperfine 
field parallel to the applied field follow. What is striking 
is that the time scales are governed by the lattice motion, 
while the hyperfine fields merely probe it. 

The temperature dependence of the homogeneous and 
inhomogeneous linewidth for Bo || c* is shown in Fig. 
both in the relaxed and completely quenched state. At 
high temperatures, the spectra are dominated by the in- 
homogeneous broadening due to the inequivalcnt Se sites. 
The line is then homogeneously broadened on the high 
temperature side of the T 2 _1 peak and is narrowed again 
on the low temperature side of the peak in the relaxed 
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state. At lower temperatures (below 10 A), the inhomo- 
geneous broadening grows significantly without evidence 
for saturation. 

The fluctuation amplitude of the hyperfine field can 
be estimated from the low temperature linewidth in the 
completely quenched state (Fig. [3J from which we ob- 
tain jSh z /2TT = O(10 A )Hz. Using this value of Sh z and 
K(6 = 90°) = 6.7 x l(r 4 , fi can be evaluated as 0.4. We 
can also estimate the characteristic time scale of the fluc- 
tuations near the peak to be t c = 0(^r X 10 _4 )s. On the 
low temperature side of the peak of T% the fluctuations 
are so slow as to be ineffective in dephasing the trans- 
verse nuclear spin components in a spin-echo experiment. 
That we should observe this trend above Tao should not 
be surprising because only on very slow cooling through 
the anion ordering transition results in a relaxed phase 
and superconductivity. 

It remains to explain the nature of the low tempera- 
ture inhomogeneous broadening in Fig. 0] which has a 
similar angular dependence as the homogeneous broad- 
ening peaking near 30A". Instead of two discrete peaks in 
the spectrum, there is a single broad feature. As a result 
of the anion ordering, there is a distinction between ad- 
jacent chains, say "A" chains and "B" chains. Zone fold- 
ing produces two bands crossing the Fermi energy, where 
the degeneracy at the zone face near (±fcp, ±7r/26, q c ) is 
lifted by a small gap, Aao- There is a corresponding 
difference in Fermi wavevector, 



The two bands are of predominantly "A" or "B" char- 
acter for these states near to the zone edge, but they 
are mixed if ^Akp < 2ir, with £ the correlation length 
of the counterion ordering measured along the stacks. 
Using Aao ~ 2meV££& and v F = 2 x 10 5 m/s gives 
Akp = O(10~ 3 )A~ 1 . Therefore, quite large values of £ 
are necessary to maintain the distinction between "A" 
and "B" chains. Perhaps sufficient order for observation 
of distinct chains could be achieved with slower cooling. 



Indeed, in magnetotransport studies^, much slower rates 
(~ lOmK /min.) were necessary to observe interference 
effects linked to the existence of the gap induced by the 
anion ordering 20 . 

If we presume the growth of the inhomogeneous broad- 
ening at low temperatures (Fig. is linked to the con- 
figuration of the counterions, then either it is a dynamic 
effect linked to motional narrowing associated with lat- 
tice defect mobility, or it is a static effect, where the na- 
ture of the defects in the counterion sublattice changes 
as the temperature is lowered. For the latter, an exam- 
ple is where the domain boundary between "A" and "B" 
stacks broadens upon cooling, perhaps as a consequence 
of thermal contraction of the lattice. Experiments sen- 
sitive to spectral diffusion could distinguish between the 
two mechanisms for the broadening. 

Finally, we make some general comments on our ob- 
servations in the context of charge ordering in organic 
conductors. In our view, the reason for the very slow 
fluctuations evident here are the large masses and mo- 
ments of inertia involved in the anion ordering transition. 
As a result of the coupling of the electronic states to the 
lattice potential, this becomes apparent in the spin-echo 
decay and inhomogeneous line broadening at low tem- 
peratures. Very similar results for spin-echo decays and 
inhomogeneous broadening were recently reported in 13 C 
NMR experiments on 6>-(BEDT-TTF) 2 RbZn(SCN)42i, a 
quasi-two dimensional organic conductor. In that case, 
there is strong evidence that the slow fluctuations are 
directly associated with a charge-ordering phase tran- 
sition, also linked to metal-insulator transition, for the 
1 /4- filled system. We suggest that the slow fluctuations 
are evidence that coupling to the lattice is an impor- 
tant component of the CO transition in that system, and 
perhaps all of the quasi-two dimensional organic systems 
undergoing charge order symmetry breaking. 
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